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In selective withdrawal, fluid is withdrawn through a nozzle suspended above the flat interface 
separating two immiscible, density-separated fluids of viscosities v^ VVRr and vio-wer = \vupper- At low 
withdrawal rates, the interface gently deforms into a hump. At a transition withdrawal rate, a spout 
of the lower fluid becomes entrained with the flow of the upper one into the nozzle. When A = 0.005, 
the spouts at the transition are very thin with features that are over an order of magnitude smaller 
than any observed in the humps. When A = 20, there is an intricate pattern of hysteresis and a spout 
appears which is qualitatively different from those seen at lower A. No corresponding qualitative 
difference is seen in the hump shapes. 

PACS numbers: 47.20.Ma, 68.05.-n, 47.55.N-, 47.20.Dr 



Fluids change shape easily in response to stress. In 
some cases, the change is so dramatic that the fluid inter- 
face changes topology. An example is fluid entrainmcnt 
by selective withdrawal. Here, a change in topology oc- 
curs when shear stresses cause the initially smooth inter- 
face between two immiscible fluids to erupt into a spout 
that pierces through the entraining fluid. If the shear is 
reduced, the spout collapses. In this paper, we address 
the nature of this topological change by characterizing 
the spout shapes near the transition. 

It is tempting to think of such topological changes as 
analogous to thermodynamic phase transitions. Some 
fluid transformations, such as drop breakup, proceed 
as some physical dimension approaches zero, becoming 
much smaller than the macroscopic dimensions of the 
flows[l|lH|. Such a separation of length scales of- 
ten leads to universal behavior as the dynamics approach 
a singularity where physical quantities diverge 0, 0, Q- 
This resembles second-order or critical behavior. While 
this framework is appealing, it is imperfect. It was re- 
cently discovered that "critical" fluid transitions can ex- 
hibit a broader set of behaviors than their thermody- 
namic counterparts, and that not all fluid singularities 
obey universal dynamics 0, 0- Other fluid transforma- 
tions, such as selective withdrawal, proceed via a dis- 
continuous jump and thus resemble first-order hysteretic 
transitions. This paper explores the nature of this type 
of topological change. 

In our experiments, schematically shown in Fig. 1, a 
nozzle is suspended a height S above the flat interface 
between two immiscible, density-separated fluids. Fluid 
is withdrawn at a flow rate Q through the nozzle, initially 
deforming the interface into an axisymmetric hump. If 
S is held fixed and Q is increased, the hump will grow 
sharper until, at a transition flow rate, a spout of the 
lower fluid becomes entrained with the flow of the upper 
fluid into the nozzle. There is hysteresis in this transition, 
and once the spout is formed, decreasing Q past a second, 
lower transition flow rate causes the spout to collapse. 

We observe different behavior depending on whether 
the selective withdrawal transition is approached from 
low Q, in the hump state, or from high Q, in the spout 
state. Previous experiments [TTL IT^ and simulations [l3j 
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FIG. 1: Experimental setup. Two immiscible fluids are den- 
sity separated in a tank 20 cm x 20 cm x 30 cm. Fluid is with- 
drawn through a nozzle a distance S above the unperturbed 
interface using a miniature gear pump. Before reaching the 
pump, the lower fluid sediments out in an airtight separa- 
tion tank, which also serves to damp pressure variations in 
the flows. The upper fluid is deposited in a reservoir and si- 
phoned back to the main tank, maintaining the upper fluid at 
a constant depth of 15-20 cm. The apparatus is illuminated 
from the rear, and a CCD camera images the interface which 
is then traced using the programs ImageJ (NIH) and Pro Fit 
(Quansoft). 



have described the approach to this transition from the 
hump state, analyzing it either as a weakly first-order 
thermodynamic phase transition^, or as a saddle- 
node bifurcation [Fj ■ 

In this paper we analyze the approach to the transition 
from the spout. We focus on two systems with different 
ratios A between the upper and lower fluid viscosities, 



and vie 



respectively. When A 



0.005, the spouts exhibit length scales near the transition 
that are an order of magnitude smaller than any observed 
in the humps. These small scales allow the entrainment 
of extremely thin fluid threads, which arc potentially use- 
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FIG. 2: Spout shapes at A = 0.005. The upper fluid is 



heavy mineral oil with v u 



= 198 cSt, p u 



= 0.87 



g/ml. The lower fluid is deionized water with vi ow<ir = 1.000 
cSt, piower = 0.998 g/ml. The surface tension between 
the fluids is 7 = 35 dyne/cm^^l- ( a ) Image of a typical 
spout. (b) z versus r, shown for three spouts with S = 0.51 
cm: (Q-Q*)/Q* = 0.019 (black circles) , (Q-Q*)/Q* = 0.65 
(gray diamonds), (Q — Q*)/Q* = 2.34 (open squares). The 
outlines are averaged for every three pixels in r and z. S 
changes by less than 0.01% while acquiring each data point. 
Inset shows 2 versus log(r). (c) z versus n ax is shown for the 
same values of S and Q as in (b). 



ful in many applications, such as the coating of biological 
tissue for the purposes of transplantation. [3 113 O n the 
other hand, when A = 20, we observe a much thicker and 
qualitatively different spout. By contrast this increase in 
the viscosity ratio does not cause a corresponding quali- 
tative change in the hump state [T3 | . In addition, at this 
value of A two distinct types of spouts appear at the same 
flow rates but in different hysterctic regimes. 

We first examine the spout shapes near the transition 
when A = 0.005. A typical spout profile is seen in Fig. 
2a. Gravity forces the interface to be horizontal far from 
the nozzle, and the entraining flows force it to be vertical 
inside the nozzle. We analyze the shape of the spatial 
region that connects these two asymptotic behaviors. As 
Q is varied, the shape of this connecting region changes, 
as shown in Fig. 2b. At low Q, the spout resembles a thin 
thread attached to a broad base, with a localized region 
of high curvature connecting the two structures. When 
Q is increased, the spouts become smoother, approaching 
a nearly logarithmic shape as shown in the inset where z 
is plotted versus logr. 
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FIG. 3: Evolution of profile parameters as flow rate is de- 
creased towards the transition, (a) Zp versus (Q — 
Q*)/Q*. (b) R„ versus (Q - Q*)/Q*. (c) r v versus (Q - 
Q*)/Q* ■ In each case, three values of S are shown: S — 0.51 
cm (black circles), S — 0.44 cm (gray diamonds), and S = 
0.32 cm (open squares). Dashed lines are power law fits. 



We characterize these profiles by calculating the axial 
and azimuthal curvatures 
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Very near the transition to the hump, the spout exhibits a 
sharp peak in K ax , as seen in Fig. 2c. As Q increases and 
the spouts become smoother, the peak in n ax disappears 
below our resolution. Wc characterize the evolution of 
the spout as the transition is approached by examining 
the following length scales: z p (the z location of the peak 
in K ax ), R p = 1/K ax (z p ), and r p = r(z p ). 

In Fig. 3, we plot these profile parameters as func- 
tions of the flow rate Q for three values of S. Near the 
transition, z p approaches a minimum value z m i n as Q 
approaches a value Q*. Using z m in and Q* as fitting pa- 
rameters, we find that plotting log (z p — z m i n )/S versus 
log (Q — Q*)/Q* collapses the data onto a single curve, 
as shown in Fig. 3a, whereas the data for R p and r p 
collapse without scaling by S, as seen in Figs. 3b and 
3c. We search for power-law scaling in the same way as 
was done in the analysis of the approach to the tran- 
sition from the hump state0. For ((Q - Q*)/Q* < 
0.6), (z p - z mm )/S oc ((Q - Q*)/Q*)°- 50±0 - 08 . Us- 
ing Q* determined from this fit, we find that R p oc 
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{(Q - Q*)/Q*)°- 45±015 . In this case, our fitting range 
is smaller: (Q - Q*)/Q* < 0.1. Although we are able to 
place an upper limit on the point closest to the transi- 
tion, our imaging does not allow us to determine R p more 
precisely. Finally, for (Q — Q*)/Q* < 1.5, r p remains con- 
stant: r p = 22 fim ± 4 /im. All three quantities show 
large departures from power-law behavior far from the 
transition. 

Previous studies of the approach to the transition from 
low Q found that the radius of curvature of the hump tip 
could be much smaller than the other lengths character- 
izing the system. For A = 0.005, the largest separation of 
length scales {h c K)h ump was 10(we estimate h c from the 
largest h max value shown) [l2^ . Applying a similar analy- 
sis to the approach from high Q, we find that z m i n /R p is 
40. We note, however, that R p is not the smallest length 
scale observed for this system. The azimuthal curvature 
r p is significantly smaller than R p , so that z m i n /r p is 
greater than 150. 

When the viscosity ratio is increased to A = 20, we 
observe a qualitatively different spout. These spouts are 
very thick, and the flows penetrate deeply into the bulk 
of the lower fluid. This is shown in Fig. 4a, where the 
flows, obtained by dyeline tracing inside the lower fluid, 
are superimposed on a typical spout image. These flows 
are distinct from those seen when A = 0.005, where the 
flows are primarily along the interface^5|. We therefore 
refer to the broad spouts as "bulk-flow" spouts. We note 
that the bulk-flow spout can be so thick that it fills the 
entire orifice when the nozzle diameter, £>, is small. We 
restrict our analysis to D large enough that this does not 
occur. 

At A = 20, a spout closely resembling those seen at 
small A also exists in a narrow hysteretic region at the 
same values of S and Q as does the bulk- flow spout. We 
are able to observe these spouts by carefully manipulat- 
ing the flow rate. At the transition flow rate, an unstable 
thread of the lower fluid becomes entrained in the upper 
one. After a time lag, which can be as long as several 
seconds, this initial spout begins to widen, culminating 
in the steady bulk-flow spout. However, if Q is rapidly 
reduced before the initial spout widens significantly, a 
stable spout as much as two orders of magnitude thin- 
ner than the bulk-flow one is produced. This thin spout 
exhibits flow patterns in the lower fluid similar to those 
seen at A = 0.005. 

The profiles of these two spouts are distinctly different, 
as seen in Fig. 4b. In the inset to Fig. 4b, the bulk- 
flow profiles are fit to z sx e~ r l TQ . The data deviates 
very slightly from this form very near the straw entry as 
well as near the horizontal interface, but is in excellent 
agreement in the intermediate region. This fit is used to 
calculate n ax {z) for these spouts. The decay length, tq, 
changes only slightly as a function of D: for D = 0.40 cm, 
r = 0.54 cm ± 0.02 cm, while for D = 0.80, r Q = 0.44 cm 
± 0.02 cm. In contrast to this exponential, the hysteretic 
surface-flow spout profiles are closer to logarithmic, and 
are similar to those seen in the inset to Fig 2b for A = 
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FIG. 4: Spout shapes at A = 20. The upper fluid is light min- 
eral oil with v U p pe r — 48 cSt, pupper = 0.85 g/nil. The lower 
fluid is glycerol with V\ ower = 980 cSt, pi ower = 1.26 g/ml. 
The surface tension between the two fluids was measured to 
be 7 = 27 ± 1 dyne/cm using the pendant drop method, (a) 
image of a typical bulk- flow spout. Flow lines were traced 
from dyelines in another image at the same S and similar Q 
and superimposed on the image shown, (b) z versus r for a 
bulk-flow spout (filled circles) at Q — 57.1 ml/s and S = 1.05 
cm. and a hysteretic surface-flow spout (open squares) at the 
same value of S but with Q = 54.7 ml/s. Inset shows ex- 
ponential fit to the bulk-flow spout, z = 1.26e~ r,/0 ' 54 . (c) z 
versus n ax for same data shown in (b). 



0.005. The hysteretic surface-flow spout also displays a 
sharper peak in n ax than the bulk-flow spout, as shown 
in Fig. 4c. 

A clearly defined transition exists between the two 
spouts. At a threshold flow rate, the stable surface- 
flow spout becomes unstable and rapidly widens into a 
bulk-flow spout with no stable intermediate structures. 
Thus, there are transitions between three distinct steady 
states. In Fig. 5, we show these transitions, including 
the relevant hystereses, on a phase diagram of the tran- 
sition flow rate Qt versus S. As seen in Fig. 5a, the 
transition from hump to bulk-flow spout can be fit to 
S = (0.29 ± 0.01)(Q t )°' 32±0 03 - Thc hysteresis in the 
hump to bulk-flow spout transition is large, especially at 
low S, and is observed to vary with D. Fig. 5b shows 
a narrow but well-defined region in which, depending on 
the initial state of the system, surface-flow spouts, bulk- 
flow spouts, and humps can all be observed. This region 
drops below experimental resolution at S — 35 ml/s. 

In conclusion, the independent variation of the viscos- 
ity ratio A and the flow rate Q dramatically changes the 
spout profiles. At A = 20 there is a broad bulk-flow 
spout which disappears for smaller A. Thin surface-flow 
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FIG. 5: Phase diagram for A = 20 showing the different hys- 
teresis regions, (a) The transitions from hump to bulk-flow 
spout (filled squares) and from bulk-flow spout to hump (open 
circles) are shown. The solid line is S = (0.29)Q? 32 . (b) A 
magnified view of the phase diagram at high Q t showing, for a 
representative data set, the surface-flow spout to hump transi- 
tion (open diamonds) and the surface-flow spout to bulk-flow 
spout transition (x's). The lines are guides to the eye. 



spouts, seen at low A, are only visible in a limited hys- 



teretic region at A = 20. Cohen|l2j observed no qualita- 
tive difference in the hump shapes near the transition for 
A between 10 -3 and 1.7. Similar analysis of our data ex- 
tends the range to A = 2001 • Thus, a qualitative change 
in the spout, not seen in the hump, is found by varying A. 
Decreasing the other control parameter, Q, while fixing 
A = 0.005, causes the spout to narrow drastically as the 
transition is approached. This produces a much greater 
separation of length scales in the spout than found in the 
hump. 

Unlike the hump interface, the spout interface is not 
bounded in the vertical direction. Thus, the spout has 
two asymptotic regimes which must be matched to each 
other by the dynamics: at large radius it is constrained 
by gravity to be horizontal and at large heights it is con- 
strained by the flows in the nozzle to be vertical. The in- 
crease in the separation of length scales as the transition 
is approached can be viewed as a decrease in the spatial 
width connecting these constraints. At high Q, the axial 
curvature, K ax , is nearly constant in z. As Q is decreased 
to the point where the spout disappears, the transition 
region becomes increasingly localized, leading to a sharp 
peak in n ax - However, the spatial transition region never 
collapses to zero and K ax is always cut off at a finite value. 
The degree to which the selective-withdrawal transition 
can approach a continuous one with a singularity is lim- 
ited by the minimum extent to which the zone matching 
the two constraints can shrink to zero. 
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